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I. INTRODUCTION 

One dimensional potential wells have bound states. They are solutions of the Schrodinger equation 
which satisfy appropriate boundary conditions. When two identical potential wells are very far apart, 
then each potential has the same energy levels, and each eigenstate is doubly degenerate. As the wells 
are brought closer together, there is communication between them, and each level is split into two. 
Similarly, if one has an array of many identical wells making a periodic potential, then one gets energy 
bands, which play, for example, a crucial role in determining the electronic properties of crystalline 
solids. To illustrate this band structure qualitatively, condensed matter physics texts usually treat the 
problem of a one dimensional periodic array of delta functions, called the Kronig Penney model. Here, 
one gets a transcendental equation for computing band edges. Another well-studied class of periodic 
potentials are the Lame potentials ^ [21 El 



where sn(x,m) is a Jacobi elliptic function 4 of real elliptic modulus parameter m (0 < m < 1). If a 
is any positive integer, it is well known that these potentials have only a band gaps and hence 2a + 1 
band edges which are analytically known (in principle). However, this is an exceptional example. The 
bottom line is that even in one dimension, there are very few solvable periodic potentials, and it would 
be desirable to have more, especially some with a richer spatial structure. 

We shall describe how to obtain new solvable periodic potentials via several different approaches. 
One way, for example is to expand our knowledge of Lame potentials to the wider class of associated 
Lame (AL) potentials [HI El 



where a, b are constants and without any loss of generality we consider a >b. A second way, for example 
is to exploit the fact that Jacobi elliptic functions are doubly periodic functions and consider PT- 
invariant complex periodic potentials [71 obtained from the Lame and the AL potentials by considering 
the anti-isospectral transformation x ^ ix + (3 (where /? is any nonzero real number) ^ . A third way 
is to consider a linear superposition of several Lame (AL) potentials which in view of the Landen 
transformations [HI ^1 and cyclic identities jLL, for Jacobi elliptic functions can be shown to be 
essentially equivalent to Lame (AL) potentials with modified parameters. Finally, we can further 
expand all these classes of solvable potentials by using the techniques of supersymmetric quantum 
mechanics to generate supersymmetric partner potentials |12j . 




(1) 




(2) 
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The outline of this paper is as fohows. In Sec. II we review some general properties of periodic 
potentials as well as the relevant formalism of supersymmetric quantum mechanics. In Sec. Ill we 
discuss the band edge energies and eigenfunctions for several Lame and AL potentials. In particular, for 
Lame potentials we derive a remarkable hitherto un-noticed relation connecting band edge energy 
eigenstates corresponding to modulus parameters m and 1 — m. We also describe how the Lame 
potential results can be vastly expanded to get solutions for the AL potentials ©• In particular, we 
show that for any integral a, b with a > b, there are a bound bands followed by a continuum band, out 
of which the top b bound bands are unusual in the sense that both band edges are of period 2K{m) 
(4:K(m)) when the integer a — 6 is even(odd), where K{m) is the complete elliptic integral of the first 
kind. In Sec. IV we obtain the supersymmetric partners of both the Lame and the AL potentials, 
thereby expanding the list of solvable periodic potentials. In Sec. V we show that using Landen 
transformations [TO] and cyclic identities obtained by us recently exact band edge eigenvalues 
and eigenfunctions can also be obtained in the case of linear superposition of Lame and AL potentials. 
In Sec. VI we show that a class of PT-invariant complex periodic potentials obtained from the Lame 
and AL potentials by the anti-isospectral transformation x ^ ix + (3, {(3 being any nonzero real 
number) are also exactly solvable problems with real band spectra. Finally, in Sec. VII we discuss 
the double sine-Gordon equation and show that in some special cases it has unusual band spectra, in 
that the majority of bound bands have both of their band edges of the same period. 

II. GENERAL PROPERTIES 

(a) Periodic Potentials: For a real potential with period L, one is seeking solutions of the Schrodinger 
equation subject to the Bloch condition 



where k denotes the crystal momentum. The spectrum shows energy bands whose edges correspond 
to kL = 0,7r, that is the wave functions at the band edges satisfy ip{x) = ziziplx + L). For periodic 
potentials, the band edge energies and wave functions are often called eigenvalues and eigenfunctions, 
and we will also use this terminology. 

A general property of eigenstates for a real potential with period L is the oscillation theorem. It 
states that band edge wave functions, when arranged in order of increasing energy Eq < Ei < E2 < 
E3 < E4^ < E^ < Eq < have periods L,2L,2L, L, L,2L,2L, ... 1. The corresponding number of 
wave function nodes in the interval L are 0,1,1,2,2,3,3,... and the energy band gaps are given by 




ikL 




(3) 
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Ai = E2 — El, A2 = E/^ — E^, A3 = Eq — E^, .... The oscillation theorem is very useful in identifying 
if all band edge eigenstates have been properly determined or if some have been missed, 
(b) Supersymmetric Quantum Mechanics: The supersymmetric partner potentials V±{x) are 
defined in terms of the superpotential W{x) by ^21 

V±{x) =W^{x)±W'{x), (4) 

where (by convention) the superpotential W{x) is related to the ground state eigenfunction of V-{x) 
by 

W(.)^-f^. (5) 
The corresponding Hamiltonians H± can be factorized as 

H_ = A+A, H+ = AA+ , (6) 

where 

so that the spectra of H± are nonnegative. It is also clear that on the full line, both H± cannot have 
zero energy modes, since the two functions ipQ^^ given by 

4^)(x) = exp(± fw{y)dy) , (8) 

cannot be simultaneously normalized. 

On the other hand, when the superpotential W{x) has period L \W{x + L) = W{x)\, and the 
eigenfunctions of both V±{x) must satisfy the Bloch condition it is easily shown that in the 
periodic case, irrespective of whether supersymmetry is broken or unbroken, the spectra of V±{x) are 
strictly identical jl3l I14j . Further, one has unbroken supersymmetry provided 

t 

Jo 



Wiy)dy = 0, (9) 



and in this case the partner Hamiltonians have identical spectra including zero modes and both ipQ^^ 
belong to the Hilbert space. As a result, unlike on the full line, for periodic potentials, whenever 
the condition © is satisfied, supersymmetry is unbroken and yet the Witten index (which counts the 
difference between the number of zero modes of V'o^'') is always zero. 

The condition Q is trivially satisfied when W{x) is an odd function of x and throughout this 
paper we shall only consider superpotentials W which are odd function of x. Further, using the known 



eigenfunctions ipi, \x) of V-{x), one can immediately write down the corresponding eigenfunctions 
ipi^\x) of V-^^{x). In particular, from eq. (jEJ it follows that the ground state of V^{x) is given by J2] 

4'"^(^) = -74— , (10) 

% \x) 

while the un-normalized excited states ipi^^ (x) are obtained from ipn ^ (x) by using the relation 



V^i-)(x) , (n>l) . (11) 



Thus by starting from an exactly solvable periodic potential V^{x), one gets a new isospectral periodic 
potential V+(x). 

A few years ago, the concept of self-isospectral periodic potentials has been defined and developed 
in detail 13 . A one dimensional potential V-{x) of period L is said to be self-isospectral if its 
supersymmetric partner potential V+(x) is just the original potential upto a discrete transformation - 
a translation by any constant amount, a reflection, or both. A common example is translation by half 
a period, in which case the condition for self-isospectrality is V+(x) = V-{x — L/2). In this sense, any 
self-isospectral potential is rather uninteresting, since the application of supersymmetry just yields a 
discrete transformation and basically nothing new. 

III. LAME AND ASSOCIATED LAME POTENTIALS 

(a) Lame Potentials: The Lame potentials as given by eq. (pQ) have a period L = 2K{m). Their 
name comes from the fact that the corresponding Schrodinger equation (with Ti = 2m = 1) 

— r -|- a{a + l)msn (x, mm = Eip , (12) 

is called Lame equation ^ I21- It is well known that for any integer value a = 1, 2, 3, . . ., the corre- 
sponding Lame potential has a bound bands followed by a continuum band 1^ and the 2a + 1 band 
edge energy eigenstates are analytically known in principle. We now obtain remarkable new relations 
(valid for any integer a) relating the band edge energy eigenvalues and eigenfunctions at two values 
m and 1 — m of the modulus parameter. 

We start from the Schrodinger equation dJ). On using the relation |^ 



"v/m sn(x, m) = — dn[zx -|- K'{m) + iK{m), 1 — m] , (13) 
and defining a new variable y = ix + K'{m) + iK{m), the Schrodinger eq. H12|) takes the form 

- V"(y) + a{a + 1)(1 - m)sn2(y, 1 - m)i;{y) = [a{a + 1) - E]iP{y) . (14) 



On comparing eqs. (\12\i and (|14|) we then have the desired relations relating the eigenstates when the 
modulus parameter is m and 1 — m: 

Ej{m) = a{a + l) - E2a~jil- m) , '4)j{x,m) ^2a-j{yA - m) , j = 0, 1, 2a . (15) 

Thus, for any integer a, at m = 1/2, one has remarkable relations 

Ej{m = 1/2) + E2a-j{m = 1/2) = a{a + 1) , Ea{m = 1/2) = a{a + l)/2 , 

ipj{x,m = 1/2) oc tp2a-j{y,m = 1/2) . (16) 

On using 

dn(2;, m) = \/l — msTi[ix + K'{m) + iK{m), 1 — m] , 

\/rncn{x, m) = — mcn[ix + K' {m) + iK{m), 1 — m] , (17) 

one can immediately and explicitly verify the relations (jlSp and ()16() for a = 1,2,3,4. 

For the a = 2 case, the Lame potential has 2 bound bands and a continuum band. The energies 
and wave functions of the five band edges are well known ^ |5] . The lowest energy band ranges from 
2 + 2m — 25 to 1 + m, the second energy band ranges from 1 + 4m- to 4 + m and the continuum starts 
at energy 2 + 2m + 25, where 6 = \/l — m + m^. For subsequent application of the super symmetric 
quantum mechanics formalism, it is convenient to have a potential V^{x) whose ground state energy 
is zero. This is easily accomplished via a re-defined potential with the ground state energy subtracted 
out. For the problem under consideration, one has V- (x) = 6msn^ {x, m) — 2 — 2m + 25. This potential 
V-{x) is plotted in Fig. 1. The eigenstates of all the band edges are given in Table 1. Note that in 
the interval 2K{m) corresponding to the period of the Lame potential, the number of nodes increases 
with energy, in agreement with the oscillation theorem. From Table 1, it is easy to check that, as 
expected, the relations (|T5|) and (fT6|) are indeed true. 

We might add here that in case a is not an integer, then for Lame potentials one has an infinite 
number of bands and band gaps for which (to the best of our knowledge) no analytic results are 
available. However, when a is half-integral, then (2a+l)/2 mid-band levels are in principle analytically 
known and each of them is doubly degenerate ^ |2] . 

(b) Associated Lame Potentials: We now expand our discussion to the band edges and wave 
functions of a much richer class of periodic potentials given by eq. ((2j called associated Lame (AL) 
potentials, since the corresponding Schrodinger equation is called the associated Lame equation 
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More precisely, we often refer to the AL potential of eq. as the {p, q) potential where p = a{a + 1) 
and q = b(b + 1) and note that (p, 0) potentials are just the ordinary Lame potentials. The AL 
potentials © can also be re-written in the alternative form V{x) = pm sn^(a:) + qm sn'^(x + K{m)) 
|1J|S|. Clearly, the potentials {p, q) and {q,p) have the same energy spectra with wave functions shifted 
by K{m). Therefore, it is sufficient to restrict our attention to p > q, i.e. a > b. 

In general, for any values of p and q, the AL potentials have a period 2K{m), but for the special 
case p = q, the period is K(m). From a physical viewpoint, if one thinks of a Lame potential {p, 0) as 
due to a one-dimensional regular array of atoms with spacing 2K{m) and "strength" p, then the AL 
potential (p, q) results from two alternating types of atoms spaced by K{m) with "strengths" p and q 
respectively. If the two types of atoms are identical [which makes p = q\, one expects a potential of 
period K{m). 

Extrema (defined for this discussion as either local or global maxima and minima) of AL potentials 
are easily found by setting dV{x)/dx = 0. Extrema occur when sn(x) = 0, or cn(x) = 0. Also, for 
fixed values of q and m, there are additional extrema if p lies in the critical range 

q{l — m) < p < q/{l — m) . 

The AL potentials for q = 2,m = 0.5 and several values of p are plotted in Fig. 2. In the critical 
range of p values 1 < p < 4 , one expects additional extrema, and these are clearly seen, 
(c) Parabolas of Solvability: The AL equation is 

2,/, ^^2/ 



+ [a{a + l)m sn\x) + b{b + l)m — - E]^ = . (18) 



dx^ dn^(x) 
On substituting ^(x) = [dn(x)]^^ y(x) , it is easily shown that y(x), satisfies the Hermite elliptic 
equation pp. On further substituting sn(x) = sint , y{x) = z{t) , one obtains Ince's equation 

{l + Acos2t)z"{t)+Bsm2tZ'{t) + (C + D cos2t)Z{t) =0, (19) 

where 

m (26 — l)m X — {a + b){a + 1 — b)m 

2 — m 2 — m 2 — m 

2-m ' ' ^ ' 

which is a well known quasi exactly solvable equation [T]. In particular if a -|- 6 -|- 1 = n(n = 1,2,3,...) 
then one obtains n solutions, which are given in Table 2. In particular, for any given choice of 
p = a{a + 1), Table 2 lists the eigenstates of the AL equation for various values of g = b{b + 1). 



For q = a(a— 1), there is just one eigenstate with energy ma^ and wave function ijj = dn°'{x). Since 
the wave function has period 2K{m) and is nodeless, this is clearly the ground state wave function of 
the (a(a + l),a(a — 1)) potential for any real choice of the parameter a. The equations p = a{a + 1) 
and q = a{a — 1) are the parametric forms of the equation of the parabola (p — g)^ = 2(p + q), which 
is plotted in Fig. 3 and labeled PI. For any point on the parabola, one knows the ground state wave 
function and energy Eq = ma^ . The parabola PI includes the points (2,0) and (6,2). 

For g = (a — l)(a — 2), we see from Table 2 that two eigenstates with energies 1 + m{a — 1)^ and 
l + ma? are known. Since they have period AK{m) and just one node in the interval L = 2K{m), they 
must correspond to the first and second band edge energies Ei and E2 of the (a(a + 1), (a — l)(a — 2)) 
potential. Eliminating a from the equations p = a{a + 1) and q = (a — 1) (a — 2) gives the "parabola 
of solvability" {p — q)^ = 8{p + q) — 12, which is plotted in Fig. 3 and labeled P2. This parabola 
includes the points (2,0) and (6,0) which correspond to Lame potentials. Similarly, the parabolas of 
solvability Pn (n = 0, 1, 2, ...) corresponding to q = (a — n + l)(a — n) in Table 2 are plotted. Note 
that n eigenstates are known for any point on the parabola of solvability Pn. 

All Lame as well as AL potentials for which a, b are unequal integers, have two parabolas of 
solvability passing through. This provides a good understanding of why these are completely solvable 
problems. For instance, the (2, 0) potential is at the intersection of parabolas PI (1 known state) 
and P2 (2 known states), thus giving the 3 known band edges. Similarly, the (6,2) AL potential, 
lies on parabola PI (1 known nodeless state of period 2K{m)) and parabola P4 (4 known states of 
period AK{m), two with 1 node and two with 3 nodes). Since we know from the oscillation theorem 
that 2 states of period 2K{m) are missing, so it would appear that this is an example of a quasi 
exactly solvable potential. However, using a well known theorem about Ince's equation pi we discuss 
below that the two "missing" states are degenerate and have no band gap between them. Hence, AL 
potentials are also exactly solvable periodic problems with a finite number of band gaps when both a 
and b are unequal positive integers. 

Other fully solvable examples correspond to AL potentials with a = b = integer. For example, the 
(2,2) potential has period K{m). It lies on parabola P3 (3 known states) and the band edge periods 
are K{m),2K(m),2K{m). We shall discuss these examples in some detail in Sec. V where we show 
that in view of the Landen transformation formulas for Jacobi elliptic functions, these problems are 
essentially related to Lame potentials with integer a and hence exactly solvable. 
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(d) Exact Results: Several exact results are known in the literature about the Ince's equation ()19() . 
In particular, it is well known that the system satisfying Ince's eq. (|19|) has at most j + 1 band gaps 
of period 7r[27r] in case the polynomial Q{iJ,)[Q*{fi)] defined by 

Qifi) = 2Afi^ - B^i- D/2, Q*{ii) = 2A{ii-l/2f - B{n- 1/2) - D/2, (21) 

where A, B, D are as given by eq. H20|) has either non-negative integral roots, the highest of which 
is J, or negative integral roots, the lowest of which is —j — 1. The fact that the Lame potentials (|T]) 
have only a band gaps when a is any positive integer is easily understood from here. We now show 
[S] that on applying this theorem, one can draw the following conclusions about the associated Lame 
potentials: 

1. a, h unequal positive integers: In this case one can show that for a > b, there are only a bound 
bands followed by a continuum band out of which the lowest a — b bands are normal bands with 
one band edge wave function of period 2K{m) and the other with period 4:K{m), while the top 
b bound bands are unusual in that both of their band edges are of period 2K{m)[4:K{m)\ in case 
a — 6 is an even [odd] integer. Further, all the 2o + 1 band edge eigenvalues and eigenfunctions 
are in principle known analytically. As an illustration, in Table 3 we have given all five band 
edge energy eigenstates of the (6, 2) potential. 

2. Both a, b being half-integral: In this case there are infinite number of bands and band gaps. 
However, but for the lowest a — b bands, the rest are unusual in that both of their band edges 
are of period 2K{m)[AK{m)] in case a — 6 is an even [odd] integer. Further, in this case a — b 
band edge energy eigenstates of period 2K {m)[4:K {m)] are analytically known in principle, in 
case a — 5 is an odd [even] integer. Besides, one also analytically knows the energy eigenvalues 
and eigenfunctions of 6 + 1/2 mid-band states of period 2K{m)[AK{m)] each of which is doubly 
degenerate in case a — 6 is an odd [even] integer [g] . 

3. Either a + b or a — b integral: When neither a nor b is integral or half-integral but either 
their sum or difference is an integer, then one again has infinite number of bands and band 
gaps. Further, one can show that if either a + b = 2N or a — b = 2N + 1, then there are at 
most -|- 1 band gaps of period 2K{m) and one analytically knows the energy eigenvalues and 
eigenfunctions of 2N + 1 band edges of period 2K{m). On the other hand, if a -|- 6 = 2A^ -|- 1 
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or a — b = 2N, then there are at most N band gaps of period 4K{m) and one has analytical 
expressions for 2N band edges of period 4K{m). 

Finally, just as a+ 1/2 mid-band states, with each being doubly degenerate, are analytically known 
in case parameter a of the Lame potential is half integral, one can show that a similar number 
of doubly degenerate mid-band states are also known for every integral value of b in the case of AL 
potentials © with half-integral a [H]. 

IV. SUPERSYMMETRIC PARTNER POTENTIALS: 

(a) Lame potentials: Let us first apply the supersymmetric quantum mechanics formalism as 
explained in Sec. II to the Lame potentials when a is positive integer. Since analytic solutions 
are known for integer values of a :2 , the supersymmetric partner potentials can be readily computed 
using eqs. (0]) and ©. We first discuss the results for small integer values of a, and then present some 
eigenstate results for arbitrary integer values of a. 

In order to use the supersymmetry formalism, we must shift the Lame potential by a constant to 
ensure that the ground state i.e. (the lower edge of the lowest band) has energy E = 0. For a = 1, one 
has V-.{x) = 2msn'^(x) — m, ipQ ^ = dn(x) and the superpotential is W = msn(3;)cn(x)/dn(x). The 
partner V-^-{x) turns out to be just V- {x — K(m)), so that this is an example of self-isospectrality. For 
a = 2, the potential is 



V^{x) = 6msn^{x) -2-2m + 26 , 5 = Vl- m + m'^ , (22) 

with a unnormalized ground state wave function -^q \x) = l-|-m-|-(5 — 3msn^(a;) The corresponding 
superpotential is 

W = —6msn{x)cn{x)dn{x)/ipQ \x) , (23) 
and hence the partner potential V+(x) for the potential V^{x) given in eq. is 

[1 + m + — 3msn^(xjJ^ 

Although supersymmetry guarantees that the potentials V± are isospectral, in this example they are 
not self-isospectral. Therefore, V+(x) as given by eq. H24|) is a new periodic potential which is strictly 
isospectral to the potential (|22jl and hence it also has 2 bound bands and a continuum band. In 
Fig. 1 we have plotted the potentials V±{x) corresponding to a = 2 for m = 0.8. Using eqs. (|10() 
and Hll|) and the known eigenstates of VL(x), we can immediately compute all the band-edge Bloch 
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wave functions for V+(x). In Table 1 we have given the energy eigenvalues and wave functions for the 
isospectral partner potentials V±{x). In summary, for integral a, Lame potentials with a > 2 are not 
self isospectral. They have distinct supersymmetric partner potentials even though both potentials 
have the same (2a + 1) band edge eigenvalues. 

(b) Associated Lame potentials: It is easily checked from Table 2 that the solution corresponding 
to q = a{a — 1) as well as one of the q = (a — 2) (a — 3) solutions are nodeless and correspond to 
the ground state. Hence, for these cases, one can obtain the superpotential and hence the partner 
potential V+. For example, consider the case of p = a{a + 1), g = a{a — 1). In this case W is given by 
W = amsn(x)cn(x)/dn(x) , so that the corresponding partner potentials are 

V± = {a:t l)am — ~^ 'ma{a =F l)sn^(x) — ma?' . (25) 
dn (x) 

These partner potentials are self-isospectral and so supersymmetry yields nothing new. As an illus- 
tration, the potential (6,2) as discussed in Table 3 is a self-isospectral potential and supersymmetry 
yields nothing new. 

Let us now consider the partner potential computed from the ground state for the p = a{a + 
l),g = (o — 2)(o — 3) case. Here ipoix) = [m[a — 1) — 1 — 5' + m(2a — l)sn^(x)](dn(x))'^~^, where 
5' = ^JY^^rrl^\^m^(l^^T^ . The corresponding superpotential W turns out to be 

^ m{a — 2)sn(x)cn(x) 2m(2a — l)sn(x)cn(x)dn(x) ^^g^ 

dn(x) [m(l — a) — 1 — 8' + m{2a — l)sn2(x)] 

Hence the corresponding partner potentials are 

V- = mo(a + l)sn2(x) + m{a - 3)(a - 2) - 2 - m{a^ -2a + 2) + 25', V+ = -V- + 2W^. (27) 

dn (x) 

These potentials are not self-isospectral. Thus one has discovered a whole new class of periodic 
potentials V+(x) for which three states are analytically known no matter what a is. In particular, 
the energy eigenfunctions of these three states are easily obtained by taking the corresponding energy 
eigenstates of V- as given in Table 2 and using eqs. and (|llj) . 

We might add here that applying second order Darboux transformations to the Lame potentials, 
Fernandez et al jl5j have obtained an interesting nonlocal effect, by which the transformed potential 
becomes an exact or approximately displaced copy of the original one. 
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V. SUPERPOSITION OF LAME AND ASSOCIATED LAME POTENTIALS 

Using the recently discovered cyclic identities 11 and Landen transformation formulas [9} llOj for 
Jacobi elliptic functions, we shall now further expand the list of analytically solvable periodic potential 
problems. In particular, we obtain the exact band edge eigenstates for potentials obtained via a certain 
specific type of linear superposition of Lame as well as AL potentials. 

Superposition of Lame potentials: Consider the following linear superposition of p Lame potentials 
with translated arguments: 

p 

m) = a(a + l)m ^ sn^(2;j, m) , (28) 
i=i 

where Xj = x + 2{j — l)K{m)/p. Note that for the special case of p = 2, this corresponds to the AL 
potential with a = b. On expressing sn^(x,m) in terms of dn^(x,m), the corresponding Schrodinger 
equation is given by 

p 

-<(x) - [a(a + l)^dn2(xj,m)]V'(x) = [E - a{a + l)p]ij{x) . (29) 



On making use of the Landen transformation formula [HI llUj 

P 1 X 

^ dn(xj, m) = — dn(— , rh) , (30) 

where 



a a 



1 " 



J2M2ij-l)Kim)/p,m), (31) 

and 



m = {m- 2)a^ + lo? ^ dn^(2(i - \)K(m) /p, m) , (32) 



one can rewrite the Schrodinger eq. (|29() in the form 

„ a{a + l)m ^ x . 

— ip (xj H 5 sn { — ,mm{x) 

a 



S-a(a + l)(p + 2Arf-^) 



V^(x) . (33) 



Here = X]j<fc=i dn(xj, m)dn(x/c, m). Using recently discovered cyclic identities, Ad can be shown 
to be HU 

p (P-2)/2 

Ad = ttVI — m + p ^ [dn{2jK{m)/p,rn) — cs{2jK{m)/p,m)Z{2jK{m)/p,m)] , p even, 

(p-l)/2 

= p ^ [dn{2jK{m)/p,m) — cs{2j K{m)/p,rn)Z {2jK{m)/p,m)] , p odd, (34) 
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where cs(x,m) = ££(££B) while Z(x,m) is the Jacobi zeta function W. On making the transformation 

\ ' sn(x,m) ' V ' / - . o 

X = ay, the Schrodinger eq. can be rewritten as a Lame equation but for the modulus parameter 
m and with energy ml: 

- + [a{a + l)msn2(y, m)]V'(y) = E^^'^ {rn)i){y) . (35) 

Here the true energy eigenvalues E{m) of the superposed potentials (|28jl are related to the eigenvalues 
E^^\m) of the Lame potentials by 

Ef\m) 1 
Ej{m)= ' : ' +a{a + l)[p + 2Ad-^], (36) 

where a and are as given by eqs. H31|) and H34() respectively while m and m are related by eq. ()32() . 
Similarly, the true eigenfunctions of the superposed potentials (|28|) are related to those of eq. (|35|) by 

■ipj{x,m) oc — , rh) . (37) 



As an illustration, for p = 2 (i.e. a = b AL case), it is well known [01 Ell that 

1 , , (1 - VI - m)^ / X 

a = -- — 7=^=, ^d = Vl-m, m = , (38) 

1 + A/l — (1 + Vl — m)"^ 

and hence the energy eigenvalues of the a = b AL potentials (|28() (with p = 2) are given in terms of 
those of Lame potential (but with the modulus parameter m) by 

a'^E^^^\m) = Ef\m) + a[a+l)^/^. (39) 

For a = 1, 2 we can immediately verify that the energy eigenvalues and eigenfunctions of the (2,2) and 
(6,6) potentials as calculated from here are identical (as they should be) to those calculated by us by 
an entirely different method where no use of either Landen transformations or the cyclic identities was 
made 5 . 

In fact we can relate the energy eigenstates of the a = b AL potentials at different modulus 
parameters by using the relationship H15|) between the corresponding Lame energy eigenstates. For 
this purpose, it is best to reexpress the relation H39|) entirely in terms of fh by noting that on using 
eq. (|551) we have 

m = . (40) 

(1 + ^/^)2 ^ ' 

Using eqs. (|38|) to (|ifll) in eq. (|T5|) we find that the AL eigenvalues are related by (where for simplicity 
we have replaced fh everywhere by m) 

(1 + v^)^ ^AL(^^)^ (1 + '^Y^t El^_^{m2) = aia+l)[l + V^+VT^] , j = 0,1, ...,2a , (41) 
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where 



4:^/m 4-^/1 — m 



For o = 1,2 we immediately verify that these relations are indeed true. Further, the corresponding 
eigenfunctions are related by 

ipj{x,mi) (xip2a-j{yA-m) , (43) 

where y = j^^^ + K'(m) + iK{m). 

Similarly, for p = 3, it is well known that 

1 

« = ij^2q ' = + " "^ (1 _^g)2(-x + 2g)2 ' 9 = dn(2K(m)/3, m) , (44) 

and hence the energy eigenvalues of (|28|) (with p = 3) are given in terms of those of Lame's with 
parameter rh by 

E{m) = (1 + 2qfE{m) + 2a{a + 1)(1 - q^) . (45) 

Thus by using Landen transformations and recently discovered cyclic identities one has discovered 
new solvable periodic potentials. On following the last section, it is easy to see that for any integer 
a > 2, corresponding to every superposed potential, supersymmetry will give us another new exactly 
solvable periodic potential. 

Superposition of Associated Lame potentials: Consider the the following superposition of AL 
potentials 

p p 
V{x, m) = a{a + l)m ^ sn^(2;j, m) + 6(6 + l)m ^ sn^(xj + K{m),m) , (46) 

where xj = x+2{j — l)K{m) /p. Now proceeding exactly as above and using the Landen transformation 
(|3U|) and the cyclic identity (|34() . one can show that the energy eigenvalues and eigenfunctions of the 
superposed AL potentials (|l6|l are related to those of AL potential (j2j (but with modulus parameter 
m) by 

E{m) = + [a{a + 1) + 6(6 + 1)] [p + 2Ad - \] , (47) 

X 

ip{x, m) oc ^{—,rh) , (48) 
a 

where a,rh and are given by eqs. ()31() . ()32() and (|35() respectively. Thus we have further enlarged 
the list of new exactly solvable periodic problems. Needless to say that using supersymmetry, one can 
further enlarge this list by obtaining partner potentials of the superposed potentials which will have 
the same band spectra. 
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As an illustration, for p = 2,3, the energy eigenvalues of the superposed potentials ()46|) are given 
below in terms of the corresponding AL potentials (with modulus parameter m changed to rfi): 

E{m) = [1 + Vl - mfE{m) + m[a{a + 1) + b{h + 1)] , (49) 

E{m) = (1 + 2qfE{m) + 2[a{a + 1) + h{b + 1)](1 - q^) . (50) 
where use has been made of eqs. (|5S)) and (|1H) . 

VI. PT INVARIANT POTENTIALS WITH REAL BAND SPECTRA 

We shall now show that the complex PT-invariant periodic potentials obtained by applying the 
anti-isospectral transformation ^S] a; — > ix + /? (where (3 is a. nonzero real number) to the Lame and 
AL potentials with a finite number of band gaps are also exactly solvable and have a finite number 
of band gaps. We give explicit expressions for the band edges of some of these potentials S . For 
the special case of the potential V{x) = — 2msn^(ix + f3, m), we analytically obtain the corresponding 
dispersion relation. 

Let us begin with the simple observation that if ^/^(x) is a solution of the Schrodinger equation for 
the real potential V{x) with energy E, then -(/'(ix + /?) is a solution of the Schrodinger equation for 
the complex potential —V{ix + /3) with energy —E. The new potential —V{ix + /3), generated by the 
anti-isospectral transformation x — > + /?, is clearly PT-symmetric and will be denoted by V^^'^\x). 

It is important to understand that the key point in obtaining these results is that unlike trigono- 
metric and most other periodic functions, the Jacobi elliptic functions are doubly periodic functions. 
This allows both V{x) and V^^'^\x) to be simultenously periodic, even though the periods are differ- 
ent. Note that the arbitrary nonzero constant /3 in the transformation x — > ix + /3 helps us to avoid 
the singularities of Jacobi elliptic functions. 

Let us first apply our approach to the Lame potentials On applying the anti-isospectral 

transformation x ^ ix + P 16,, where /? is any non-zero real number, it is easily shown that the 
band-edge eigenvalues and eigenfunctions of the PT-invariant potentials 

V^^'^\x) = -a{a + l)msn'^{ix + P,m) , a =1,2,3,..., (51) 

are related to those of the Lame potentials Q by [Hj 

E^^'^\m) = -E2a-j{m) , 'ip'j^'^\x,m) (x ijj2a-j{ix + P,m) , j = 0,1, ...,2a . (52) 
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It may be noted here that the PT-invariant complex potential ()51() is a periodic potential with period 
2K' {m){= 2K{1 — m)). Thus we have shown that the PT-invariant periodic potential (|51j) too has 
precisely a band gaps and hence 2a + 1 band edges at energies given by relation (|52|) . Special mention 
may be made of the remarkable fact that for any integer a, all bands and band gaps exchange their 
role as one goes from a Lame potential to its PT-transformed version. 

For any band structure problem, one important quantity is the discriminant A 1 which gives 
information about the number of band gaps as well as their widths. The question is if we can relate 
A*^^"^) with the discriminant A for the corresponding Lame potential. We now show that using eqs. 
((T3|) and this is indeed possible. In particular, on using these two equations we deduce that 

Ef'^\m) = Ej{l -m) - a{a + l) , j = 0,1, ...,2a , (53) 

and hence the corresponding discriminants are related by |H] 

A(^^)(S,m) = A[S + o(a + 1),1 -m] . (54) 

As an illustration, in Fig. 4 we plot the real and imaginary parts of the PT-invariant complex 
potential — 6msn^(ix + j3,m). Using the well known results for the Lame potential with a = 2 (see 
Table 2) and eq. the ground state (lowest band edge) eigenvalue and eigenfunction is easily 

shown to be 

■il)g{x) = sii{ix + I3,m)[l + m — 5 — 2>msv?{ix + I3,m)\ , Eg = —2 — 2m — 25 , (55) 



where 6 = \/l — m + m?. In Table 4 we have given all the five band edge eigenvalues and eigenfunctions 
of this PT-invariant potential, where we have subtracted off the ground state energy from the potential 
so that the lowest band edge by construction is at zero energy. Observe from the table that the band 
edges are both periodic as well as anti-periodic with periods of 2K'(m) and 4K'(m) respectively 
thereby showing that contrary to the suggestion of Bender et al. ^7|, the absence of anti-periodic 
band edges is not a general property of PT-invariant periodic potentials. 

Additional analytically solvable finite band gap potentials can be obtained from here by using 
supersymmetry. In particular, consider the potential v!i^'^\x) = —6msn'^{ix-\-(3,m) — Eg. Its ground 
state eigenfunction is given in eq. H55|) . and we find that the corresponding superpotential is 

TirCPTi/ N .cn(ix + j3,m)dn(ix + B,m) ^. cn(ix + P , 'm)sn(ix + P , 'm)dn(ix + /3 , m) 

'[x) = —i h oim ^- -j . (56) 

sn(ix + p, m) \\ + m — o — 3msn^(ia; -|- p, m)\ 
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Thus the supersymmetric partner potential is 

= [W^P^^Wx) + [W^P^^Ux) , (57) 

and it has the same energy spectrum as V_ (x) [see Table 4]. In this way, one has discovered yet 
another PT-invariant complex potential with a finite number of band gaps. The corresponding band 

( PT) 

edge energy eigenfunctions are immediately obtained from those of V_ by using supersymmetry 
transformations given in Sec. II. 

Yet another complex PT-invariant potential having the same band edge eigenvalues as above is 
obtained by starting from the partner potential V+{x) of the a = 2 Lame potential [eqs. H24|) and ()22() ] 
and then obtaining the corresponding PT-invariant potential by the anti-isospectral transformation 
X ^ ix + (3. Thus one has found three distinct complex PT-invariant periodic potentials having the 
same band spectra. One can show that this is true for any integral a > 2. On the other hand, for 
a = 1 one can show 8 that the three potentials are not distinct but are self-isospectral (i.e. differ 
from each other by a translation by a constant) jl3l I15j . 

For the special case of the a = 1 Lame potential, the dispersion relation is analytically known 
We now show 8_ that it is possible to derive the dispersion relation for the PT-invariant potential H51|) 
with a = 1. To that end we start from the Schrodinger equation: 

- ip'^x) + [l + m- 2msn'^{ix + /3, m)]'i/j{x) = E'ijj{x) , (58) 

where we have subtracted the ground state energy from the potential so that the new potential has 
zero ground state energy. On substituting y = ix + (3, eq. (|58() takes the form 

- '4'"{y) + [2msn2(y, m) - m]'4){y) = (1 - E)i){y) . (59) 

Now it is well known that two independent solutions of this equation are given by |H] 

H{ix + (3±ai)eyi.^[^{ix + (3)Z{ai)] 

= ' ^ ^ 

where H, 6, Z are Jacobi eta, theta and zeta functions while ai is related to E of eq. by 

E = ■msn'^{ai,m) . (61) 

On using the fact that while 6{ix + (3) is periodic function with period 2K'{m), H{ix + (3) is only 
quasi-periodic |0, i.e. 

H{i[x + 2K'{m)] +/3) = H{ix + /?) exp[-7r K'{m)/K{m)] , (62) 
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and using the Bloch condition, it is easily shown that for the PT-invariant complex potential © with 
a = 1, the dispersion relation is given by 

vr vr 
k = T—TTT, — z ^ iZ{ai) + i—— — -, (63) 
^2K'{m) ^ ' 2K{m) ' ^ ^ 

where ai is given by eq. 

We now turn to the associated Lame potentials given by eq. where without any loss of generality 
we consider a > h with both being positive integers. On using the anti-isospectral transformation, it 
is easy to see that the band edges of the potential © and its PT-transformed one as given by 



V{x) = -a{a + l)msn^{ix + /3, m) - 6(6 + ^ .__T^' 1\ ' (64) 



cn^(ia; + m) 
dn^(ix + /9, m) 



are again related by the relation H52|) . 

For the special case a = b, using relations (|4n|) . (|41|) and (|52j) we obtain an interesting relation 
betwen the eigenvalues of the AL potential and the corresponding complex PT-invariant potential 
given by 

E) >{mi) = { ) Ej{m2) -Aa{a + l) , (65) 

1 + ^Jm (1 + yjmy 

where tti-i 2 are as defined by eq. ()42() . Hence we can immediately relate the discriminant A for the 

a = b AL potential and the corresponding complex PT-invariant potential, i.e. we have 



A(-)(^,m,) = A[(/ + ^ fjE + ^a{a + l)il + V^+V^ 
1 + yjl-m (1 + yJmY 



VII. UNUSUAL BAND STRUCTURE FOR THE DOUBLE SINE-GORDON EQUA- 
TION 

Finally, we would like to consider the band spectra of the double sine-Gordon equation characterized 
by the periodic potential 

V{x) = 6^ sin^ 2x + 2a6 cos 2x , (67) 

where a, 6 are real, and without loss of generality a is taken to be nonnegative. This equation arises in 
several areas of condensed matter physics. It is well known that this system has an infinite number of 
bands and band gaps. Further, it is also well known that if a is an integer then a band edges of period 
7r(27r) are analytically known depending on whether o is an odd(even) integer jl81ll9j . However, what 
is not so well known (1| is that in case a is an integer, then the potential has unusual band spectra. 
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We start from the Schrodinger equation for the potential H67() . On using the ansatz 



ip{x) = exp(— 6cos 2x/2)</)(x) 



(68) 



in the Schrodinger equation, it is easily shown that (p satisfies Ince's eq. (fTU)) 



cp" (x) + 2bsm2x cp' {x) + [E - 2b{a - I) cos 2x](j){x) =0. 



(69) 



As discussed before, this is a quasi-exactly solvable equation, i.e. a band edges of period 7r(27r) are 
analytically known when a is an odd (even) integer 119j . We can now compute the polynomials 
Q{fj,) and Q*{n) defined in eq. (PT|) . and from them we conclude that if a is an odd (even) integer 
then there are at most ^^{^^) band gaps of period 7r(27r). 

It is worth pointing out here that similar conclusions can also be drawn about the complex PT- 
invariant periodic potential 



Note that this potential is invariant under the parity transformation x — > x + tt/2 followed by T. As 
has been shown in ref. when a is an even integer then PT symmetry is spontaneously broken 

and there are no quasi-exactly solvable states of period 2tt since in this case all the eigenvalues are 
now complex conjugate pairs. However, if a is an odd integer then there are a quasi-exactly solvable 
states of period vr and in this case PT symmetry is not spontaneously broken. On completing the 
above analysis, it follows that if a is an odd integer then there are at most (a + l)/2 band gaps of 
period vr. Thus in this case the band spectrum is rather unusual in that the majority of the infinite 
number of band edges are anti-periodic and have period 2tt, thereby again showing that the absence 
of anti-periodic band edges is not a characteristic feature of PT-invariant potentials. 
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Table 1: The eigenvalues and eigenfunctions for the 5 band edges corresponding to the a = 2 
Lame potential V- as given by eq. (|22|) which gives (p, = (6,0) and its SUSY partner V+. Here 
B = 1 + m + (5 and 5 = \/l — m + m?. The potentials V± have period 2K{m) and their analytic forms 
are given by eqs. ((22) and 1)241) respectively. The periods of various eigenfunctions and the number of 
nodes in the interval 2K{m) are tabulated. 



E [B - 3m sn^ {x)]iIj^^^ Period Nodes 

m + 1 + 5 -Smsn^ix) 1 2K 

26 — 1 — m cn(x)dn(x) sn{x)[6m — [m + 1)B -\- msn'^{x){2B — 3 — 3m)] 4K 1 

26-l + 2m sn(x)dn(x) cn{x)[B + msn'^{x){3 - 2B)] AK 1 

26 + 2 -m sn(2;)cn(x) dn{x)[B + sn'^{x){3m - 2B)] 2K 2 

4:6 m + 1 — 6 — 3171811^ (x) sn(x)cn(x)dn(a;) 2K 2 



Table 2: Eigenvalues and eigenfunctions for various associated Lame potentials {p, q) with p = 
a{a + 1) and q = [a — n + l)(a — n) for n = 1,2,3, .... The periods of various eigenfunctions and 
the number of nodes in the interval 2K{m) are tabulated. Here 64 = y^l — m + m^{a — 1)'^ , 6^ = 
\/4 — 7m + 2ma + m^{a — 2)^ , = ^J4 — m — 2ma + m'^{a — 1)^ , 67 = \/9 — 9m + m^(a — 2)^ . 



q E dn "(x)^/' Period Nodes 

a{a — 1) ma^ 1 

(a-l)(a-2) l + m(a-l)2 ^ 

(a-l)(a-2) l + ma2 ^ 

(a-2)(a-3) 2 + m(a2 - 2a + 2) ± 2^4 [m(2a-i)sn^(xM+m-ma±^4] 

(a-2)(a-3) 4 + m(a - 1)^ ^^ggl^ 

(a-3)(a-4) 5 + m(a2 - 4a + 5) ± 2,55 cn(x)[m(2a^l)sn2jx)-2+2>n-ma±^5] 

(a-3)(a-4) 5 + m(a2 - 2a + 2) ± 2^6 M-)lm{2a-l)snHx)-2+m-ma±5e] 



dn'*(x) 



/ ,\/ r\ in I ^2 /I I t:^ I or sn(x)cn(x-)[m(2a--l)sn2(x)— 3+2m— maztiSrl 

(g — 4)(a — 5) 10 + m(a^ — 4a + 5) =fc 207 yn \ j^^-j^-^) 
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2/r 





4/s: 


1 


AK 


1 


2K 


2,0 


2K 


2 


AK 


3,1 


4K 


3,1 


2K 


4,2 



Table 3: The five eigenvalues and eigenfunctions for the self-isospectral associated Lame potential 
corresponding to a = 2, 6 = 1 which gives {p,q) = (6,2). The potential is V-{x) = 6msn^(x) + 
2m ™a|i^| — 4m, and has period 2K{m). The number of nodes in the interval 2K{m) is tabulated. 



E ijA ^ Period Nodes 








Mx) [3"^sn2 


dn^(x) 




2K 





5 — 3m — 2\/4 — 3m 


sn(a;) 
dn(x) 
sn(a::) 
dn(x) 


(x) - 2 - V4 - 


- 3m] 




1 


5 — 2m — 2\/ 4 — 5m + 


[3msn^(x) — 


2 - m - \/4 - 


5m + m?] 


AK 


1 


5 - 2m + 2\/4 - 5m + m2 


[3msn^(x) — 

cn(a;) rr, 2 

dIIM[3"isn 


2 - m + \/4 - 


5m + m?] 


AK 


3 


5 - 3m + 2^4 - 3m 




(x) - 2 + V4 - 


- 3m] 


AK 


3 



Table 4: The eigenvalues and eigenfunctions for the 5 band edges corresponding to the PT-invariant 
potential vd^"^^ = — 6msn^(ix + /3, m) — Eg where Eg = —2 — 2m — 2\/l — m + m^ [eq. (|55j) ] . The 
potential has period 2K'{m). The periods of various eigenfunctions are also tabulated. 



E 






Period 





1 + m — (5 — 


3msn^(zx + /3, m) 


2K'{m) 


m - 2 + 2(5 


cn(ix + f3, 


m)sn(ix + /3, m) 


AK'{m) 


1 - 2m + 25 


sn{ix + /?, 


m)dn(ix + /?, m) 


AK'{m) 


1 + m + 2(5 


cn(ix + P, 


m)dn(zx + /3, m) 


2K'{m) 


4(5 


l + m + 6- 


3msn^(ix + /?, m) 


2K'{m) 
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2 8 

Figure 1: The (6,0) Lame potential V^{x) = 6msn^(x,m) — 2 — 2m + 2\/l — m + rn^ [thick hne] and 
its supersymmetric partner potential V+(x) [thin line] as given by eq. H24|). The curves are plotted for 
the choice m = 0.8 . Both V-{x) and V+(x) have their lowest band edge at zero energy. 
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5 10 15 20 

P 



Figure 3: Parabolas of solvability. This figure illustrates all associated Lame potentials {p, q) which 
are quasi solvable. Each parabola corresponds to a choice of q in Table 2. Parabola Pn is for 
q = {a — n+ l)(a — n) for n = 1, 2, 3, .... One knows n eigenstates for any point on parabola Pn from 
Table 2. 
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Figure 4: A plot of the real [solid line] and imaginary [dashed line] parts of the PT-invariant complex 
potential (|51|) for a = 2. More explicitly, the potential is — 6msn^ {ix + /3, m) + 2 + 2m + 2^/1 — m + m? 
and has its lowest band edge at zero energy. The plot is made for the choice m = 0.8 and [5 = 0.4. 
The potential has a period 2/C'(0.8) = 3.3192. 



27 



